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This talk provides a brief summary of the status of lattice QCD calculations of the 
light quark masses and the kaon bag parameter B^. Precise estimates of these 
four fundamental parameters of the standard model, i.e., m u , m^, m s and the CP 
violating parameter n, help constrain grand unified models and could provide a 
window to new physics. 



1. Introduction 

It is a great pleasure to be here to help celebrate Pran Nath's 65" 1 birth- 
day. I first met Pran and Dick in 1982 when I came to Northeastern as 
a post-doc in their group. They had just developed the first supergrav- 
ity grand unified model in collaboration with Ali Chamseddine and were 
very keen to know if grand unified theories could be formulated on the 
lattice. The first hurdle in this quest was whether supersymmetric non- 
abelian gauge theories could be formulated on the lattice. Needless to say 
we did not make much progress on that front. Over the last 22 years there 
has been significant progress on formulating QCD with chiral symmetry on 
the lattice (staggered, twisted mass Wilson 1 , overlap and Domain Wall 
(DW) fermions 2 ) and some on chiral gauge theories. However, the original 
problem of formulating supersymmetric non-abelian gauge theories on the 
lattice and exploring their strongly coupled sector using lattice simulations 
still remains. 

My three years stay at Northeastern was very productive. Pran and 
Dick gave me total freedom to pursue lattice QCD and statistical mechan- 
ics and were extremely supportive of my work. During this time I continued 
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my collaboration with Apoorva Patel at Caltech and developed a number 
of new ones including those with Steve Sharpe and Greg Kilcup at Har- 
vard, Bclcn Gavela and Rich Brower who were visiting Harvard, Gerry 
Guralnik and Chuck Zemach at Los Alamos, R. Shankar at Yale, and Bob 
Cordery and Mark Novotny at Northeastern. Being around Ali, Dick and 
Pran (in those days they seemed joined at the hip) was very inspirational. 
They would come in very early in the morning, get the coffee going, lock 
themselves in the conference room and work through the day taking only 
food and toilet breaks. Unfortunately, I could not stomach coffee and thus 
missed induction into the super world. 

In this talk I would like to present a status report on two different 
calculations that were initiated during my stay at Northeastern. The first 
is on estimates of light quark masses, which is equivalent to validating QCD 
by reproducing the hadron spectrum. The second is the calculation of the 
kaon bag parameter Bk which is illustrative of QCD correction to weak 
matrix elements. Together these calculations address fixing the values of 

four fundamental parameters of the standard model the quark masses 

m u , md, and m s and the CP violating parameter r\ in the CKM matrix. 
Since any grand unified model will have to provide information on their 
origin and values, and any extension to the standard model will impact their 
values, determining their precise values is an important step in looking for 
new physics 3 . There is therefore a deep connection between lattice QCD 
and the super world. 

This talk will focus on providing the current best estimates. For a 
background on the evolution of these calculations I refer to the proceedings 
of the yearly lattice QCD conferences 4 . Also, I will use two previous 
reviews 5,6 as starting points and update the results here. 

2. Light Quark Masses 

A self-consistent calculation of light quark masses is equivalent to validating 
QCD by demonstrating that it reproduces the hadron spectrum. The cen- 
tral question is do there exist values for the coupling constant a s {Mz) 

and the five quark masses m u , rrid, m s , m c and rrib such that lattice QCD 
reproduces the masses and decays of all hadrons? Over the last two decades 
the answer is slowly but surely converging towards YES. 

There are three methods that have been used to estimate the masses of 
the three light quarks — up, down and strange. These are chiral perturba- 
tion theory (xPT), QCD sum rules, and lattice QCD. It is instructive to 
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assess the strenghts and weaknesses of each of these methods. 

%PT is an effective theory of pseudoscalar mesons. Terms in the chiral 
Lagrangian have the same symmetries as QCD and are classified in powers 
of m 2 and p 2 . The expansion is characterized by the Gasser-Leutwyler (GL) 
coefficients 7 . Using the xPT Lagrangian one derives expressions for the 
masses and decays as an expansion in the quark masses and momenta, GL 
parameters, and an unknown scale A x pt- An example of such a relation, 
which is relevant for the discussion of quark masses, is the expansion for 
M\ tailored to staggered fermions on the lattice 8 



(M 1 K J oop ) 2 
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where the Li are the GL constants, l(M 2 ) = ln(M 2 /A 2 , PT ) are the chiral 
logarithms, and terms proportional to a 2 contain the leading discretization 
errors. 

Both xPT and lattice QCD rely on relations like Eq. 1 to relate hadron 
masses to quark masses. The advantage of lattice QCD over x?T is that 
one can carry out simulations dialing the quark masses and thus validate 
Eq. 1 over a range of quark masses. In xPT one is limited by the number 
of physical meson masses and, unfortunately, even at NLO there are more 
unknown GL coefficients than there are pseudoscalar mesons. The key point 
made by Sharpe and Shoresh 9 is that as long as one simulates with 3 flavors 
of dynamical quarks (with equal or unequal valence and sea quark masses) 
that are small enough for NLO xPT to be reliable, the GL coefficients are 
the same as in physical QCD so extrapolations to the physical point can 
be made using Eq. 1 and the GL coefficients so extracted are those of the 
physical world. 

The second limitation of xPT is that since it is based on the symmetries 
of QCD it docs not provide an absolute scale for quark masses but does 
well in predicting ratios. At the 1-loop level it provides two ratios (one, due 
to the Kaplan-Manohar symmetry, requires some additional but reasonable 
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assumptions) 



10 



2m, 



'S 



24.4(1.5) 



m u + m d 
m u 



0.553(43) . 



(2) 



m d 



Estimates of these ratios from chiral perturbation theory are more accurate 
than present lattice results. The main uncertainties in lattice calculations 
of m u and nid are due to lack of complete control over chiral extrapolations 
down to a few MeV and due to ignoring electromagnetic effects in the 
simulations. So this talk will focus on lattice estimates of m s and knowing 
it m u and rrid can be estimated using the ratios given in Eq. 1. 

Lattice simulations until 2000 were done mostly in the quenched ap- 
proximation due to limitations of computer power. In this approximation 
the effects of virtual quark loops on background gauge configurations are 
neglected. This approximation, in principle, is drastic as the quenched the- 
ory is non-unitary and one relied on phcnomenological arguments to assume 
that estimates are reliable to within 10%. Nevertheless, since calculations 
of observables on an ensemble of quenched background gauge configurations 
are, in most cases, the same as in the full theory, the community obtained 
valuable understanding and control over statistical and the following sys- 
tematic errors 

• Finite volume corrections 

• Extrapolation to the continuum limit 

• Chiral extrapolation to physical m u and m d . 

• The calculation of renormalization costants using improved pertur- 
bation theory and/or non-perturbative methods. 

What was missing was control over quenching errors and estimates of light 
quark masses using different states to set the lattice scale and fix their 
values varied by 10 — 30% 5 . 

In the last year two collaborations, HPQCD-MILC-UKQCD 11 and CP- 
PACS/JLQCD 12 , have reported results based on simulations with 2 + 1 
flavors of dynamical quarks. These simulations explore a range of quark 
mass values (as low as m s /8 in the MILC simulations) and yield 



m s (MS, 2 GeV) = 76(0) (3) (6) (0) MeV (HPQCD - MILC - UKQCD) 
m s (MS, 2 GeV) = 80.4(1.9) MeV (CP - PACS/JLQCD) (M K ) 
m s (MS, 2 GeV) = 89.3(2.9) MeV (CP - PACS/JLQCD) (M ) . (3) 
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The CP-PACS/JLQCD collaboration use the estimate from Mk for their 
central value and increase the upper limit of the error to accommodate the 
estimate from M^. Last year, based on preliminary results from these two 
collaborations, I had concluded that m s (MS, 2 GeV) = 75(15) 5 . Both 
groups have tightened control over some of their uncertainties and the new 
estimate, averaging the two results based on Mr-, is 



i.e. mostly the change is in a reduction in the error estimate. 

This value is significantly lower than the estimate from QCD sum rules 
which until 1996 was m s (~MS,2 GeV) = 125(40) MeV 10 . Since then esti- 
mates from QCD sum rules have been getting lower and tracking those from 
lattice QCD 5 . The two main uncertainties in QCD sum rule analysis are 
the quality of experimental information on spectral functions in the scalar 
and pseudoscalar channels on one side and the convergence of perturbation 
theory on the other 5 . The most hopeful channel for precise determination 
of m s is hadronic r decays. New precision data in this channel have been 
reported by the CLEO 13 and OPAL 14 collaborations. These have lead to 
better understanding of the SU(3) breaking effects in the hadronic r-decay 
sum rules and better resolution of the scalar and vector spectral functions. 
Incorporating these results Gamiz et al. derive the estimate 15 



The analysis remains sensitive to resolving SU(3) breaking, i.e. the cancel- 
lations between ud and us parts, as pointed out by Maltman 16 . Hopefully 
the precision of sum rule analysis will improve as more experimental data 
is collected. 

The bottom line is that lattice QCD has revised our thinking regarding 
light quark masses, i.e., m s is significantly lower than the estimate used by 
phenomcnologists until 1996 m s (MS,2 GeV) = 125(40) MeV 10 . A lower 
estimate has two important consequences. One, it increases the standard 
model estimate for e'/e 17 and, second, it poses a challenge to grand unified 
model builders 3 . 

2.1. Is m u = 0? 

As the quality of lattice data improve with respect to both the number of 
quark masses and lattice scales simulated one can make increasingly precise 
fits to relations like Eq. 1 . Through these fits one can extract various GL low 



m s (MS,2 GeV) = 78(10) 



(4) 



m 8 (MS,2 GeV) = 81(22). 



(5) 
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energy couplings in the chiral Lagrangian. Here I summarize the extraction 
of 2Lg — L5 by the MILC collaboration 11 who find 

2L 8 -L 5 = -0.2(1)(2) x 10- 3 , (6) 

which is significantly different from the range 

-3.4 x 1CT 3 < 2L 8 - £ 5 < -1.8 x 1CT 3 (7) 

allowed by x?T for m u — 0. So current lattice data rule out m u = which 
would have provided a convenient solution to the strong CP problem. 



3. B K 

The kaon bag parameter B K measures the QCD corrections to the weak 
mixing between K° and K°. It is defined as the dimensionless ratio of the 
matrix element of the AS — 2 effective weak operator between a K° and 
K° to its vacuum saturation approximation 

o _ (K°\s^(l- l5 )d - l5 )d\K°) 

K l(KZ\- S1 ^d\Q)(0\- S1 ^d\K°) • 1 1 

Its measurement gives a constraint in the form of a hyperbola in the f) = 
77(1 - A 2 /2), p = p(l - A 2 /2) plane 17 

fj[(l - p)A 2 r l2 So(x t ) + P (e)]A 2 B K = 0.226 

where 77, p and A are parameters in the CKM matrix. 

Over time different approaches, including chiral perturbation theory, the 
large N c expansion, QCD sum rules and lattice QCD, have been used to es- 
timate Bk ■ Current phenomenology uses the lattice result obtained by the 
JLQCD collaboration with unimproved staggered quarks in the quenched 
approximation 18 

B K (MS, NDR, 2 GcV) = 0.63(4) (9) 

or the corresponding renormalization group invariant quantity 6 

B~k = 0.86(6) (14), (10) 

where the second error is an estimate of the quenching and SU(3) breaking 
(nid and m s are degenerate in these calculations) uncertainty. 

The JLQCD calculation 18 showed that there are significant (i) a 2 dis- 
cretization errors in unimproved staggered fermions and (ii) the unknown 
0(a 2 ) corrections to the one-loop perturbative renormalization constants 
could be as large as ~ 10%. The continuum extrapolation involved a subtle 
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Figure 1. Summary of quenched estimates for Bk(M S , N DR,2 GeV) from different 
collaborations discussed in the text. The band 0.58 ± 0.04 captures the spread in esti- 
mates from recent simulations using improved fermion formulations (Domain wall (DW), 
overlap, and staggered(Stag)). 



interplay between these two effects. For this reason there has been concern 
whether the continuum extrapolation leading to Eq. 9 is under control. To 
test this, and because the calculation of Bx provides a laboratory for eval- 
uating the efficacy of improved fermion formulations, there are a number 
of new simulations with different lattice formulations. 

I will focus on three sets of new calculations, all within the quenched 
approximation, to elucidate our current understanding of the continuum 
limit, (i) Fermion formulations that incorporate chiral symmetry a la 
Ginsparg- Wilson (Overlap and Domain Wall fcrmions); (ii) improved stag- 
gered fermions; and (iii) twisted mass lattice QCD. These results are sum- 
marized in Fig. 1. 
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• The domain wall simulations by CP-PACS 19 and RBC 20 collab- 
orations agree with each other and give Bk (MS, NDR, 2 GeV) = 
0.57(3) after extrapolation to the continuum limit. 

• There are two calculations using overlap fermions: DeGrand 21 
finds B K (MS, NDR, 2 GeV) = 0.55(3) at j3 = 5.9 and 6.1 whereas 
Garron et al. 22 find B K (MS, NDR, 2 GeV) = 0.63(6) at = 6.0. 

• Improved staggered fermion simulations by Lee et al. 23 at [3 = 6.0 
yield B K (MS,NDR,2 GeV) = 0.58(2)(4). This calculation shows 
that both discretization errors and perturbative corrections to 
renormalization constants are small with HYP smeared staggered 
fermions, alleviating the two most serious systematic errors in this 
approach, and making it computationally attractive. 

• Estimate from simulations of Twisted mass QCD by Dimopoulos 
et al 24 is Bk (MS, NDR, 2 GeV) = 0.592(16) in the continuum 
limit. 

All these new estimates are, within errors, consistent and covered by 
the range 

B K (MS, NDR, 2 GeV) = 0.58(4) (11) 

shown in Fig. 1. Even though no single collaboration has obtained suffi- 
cient data to do a reliable continuum extrapolation, what has become clear 
from these calculations is that by improving the lattice formulation (and 
in some cases supplementing it with non-perturbative determination of the 
renormalization constants) the dependence on a has been reduced very sig- 
nificantly. What is less clear is whether the difference from the JLQCD 
result 0.63(4) is significant. With hindsight one can look at the data in 
Fig. 1 and conclude that the extrapolation of the unimproved staggered 
estimates to the continuum limit is at fault because the fit parameters were 
not well determined by the data. My take on this issue is that, given the 
data, JLQCD did the best possible extrapolation incorporating the lead- 
ing two corrections, a 2 and a 2 (since they used one-loop matching between 
the lattice and MS schemes), and requiring that gauge- invariant and non- 
invariant operators give the same estimate in the continuum limit. Their 
analysis including only a 2 errors gave 0.598(5). On adopting a better mo- 
tivated procedure their error estimate increased considerably and explains 
the difference at roughly ler level if indeed the final quenched value settles 
at 0.58. 

Computationally, the simplest of these approaches to extend to dynami- 
cal quarks is the improved staggered, however, there is a caveat. Dynamical 
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staggered simulations require taking the fourth root of the determinant in 
order to simulate one flavor of quarks. It has not been shown rigorously 
whether the action of the resulting theory is local or in the same universal- 
ity class as QCD. Simulations by the MILC collaboration 11 of a number of 
observables suggest that the continuum limit is approached smoothly and 
this fourth-root trick does reproduce QCD. Work on clarifying this issue is 
in progress. 

Dynamical simulations are just beginning. The Rikcn-Brookhaven- 
Columbia collaboration have presented first results based on three en- 
sembles of n/ = 2 dynamical lattices with domain wall fermions at scale 
a -1 « 1.7 GeV and quark masses in the range m s — 0.5m s 25 . Their result 
is 

B K [MS, NDR, 2 GeV) = 0.509(18) (12) 

for degenerate quarks (mj = m s ) and 

Bk (MS, NDR, 2 GeV) - 0.495(18) (13) 

with nid — 0.5m s . Compared to their quenched estimate 0.57(3) they find 
a ~ 15% decrease. It remains to be seen how much more this estimate will 
change when a third flavor is added to the simulations and ma is varied 
over a larger range. 

4. Conclusions 

Much has changed since I was at Northeastern. The campus has been 
improved beyond recognition. Many things are the same. All my friends 
and colleagues are still thriving and Pran is as productive, dedicated and 
driven as ever. He continues to be an inspiration for all. 

Lattice QCD has progressed tremendously. Gone are the days when one 
stayed up nights trying to harness all possible VAX computers running at a 
fraction of a mcgaflop to generate one quenched background configuration 
in a month. One now talks of sustained teraflops on dedicated massively 
parallel computers. As a result of this increase in computing power and bet- 
ter algorithms and theoretical understanding we are now simulating QCD 
without any approximations (with 2 + 1 dynamical flavors). Thus, from now 
on the community will provide increasingly precise estimates and hopefully 
one day soon the effort will, without doubt, validate QCD and yield a 
glimmer of new physics. 
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